Infrared Behavior of Interacting Bosons at Zero Temperature 
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We exploit the symmetries associated with the stability of the superfluid phase to solve the long- 
standing problem of interacting bosons in the presence of a condensate at zero temperature. Im- 
plementation of these symmetries poses strong conditions on the renormalizations that heal the 
singularities of perturbation theory. The renormalized theory gives: For d > 3 the Bogoliubov 
quasiparticles as an exact result; for 1 < d < 3 a nontrivial solution with the exact exponent for the 
singular longitudinal correlation function, with phonons again as low-lying excitations. 
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The problem of understanding the low-lying excita- 
tions from the ground state of an interacting Bose system 
has been one of the major problems of condensed mat- 
ter theory in the fifties and sixties. The first solution 
to this problem was given by Bogoliubov as a general- 
ized Hartree-Fock approximation [1]. Numerous papers 
were then devoted to analyze the corrections to this ap- 
proximate solution [2-12]. Apart from approximations 
showing a gap in the excitation spectrum, all attempts to 
improve the Bogoliubov approximation encountered the 
difficulty of a singular perturbation theory (PT) plagued 
by infrared (IR) divergences, due to the presence of the 
Bose-Einstein condensate and the Goldstone mode [2,3]. 

A systematic study of these IR divergences has been 
long delayed because they appeared only at intermediate 
steps of the calculations while physical quantities turned 
out to be finite [2] . That problems could arise in PT was 
originally recognised by Gavoret and Nozieres [3] ; indeed 
later [5] it was found that the exact anomalous self-energy 
(at zero external momentum) has to vanish [13], in con- 
trast with the finite value obtained by the Bogoliubov 
approximation. This result questions the validity of PT 
and requires a properly renormalized theory. 

To take care of the IR divergences at arbitrary spa- 
tial dimension d greater than one, we exploit the 
renormalization-group (RG) approach. In its standard 
application the RG approach sums up the singularities 
of PT and provides the power-law behavior of physical 
quantities which is characteristic of critical phenomena. 
Here we deal instead with a stable superfluid phase, for 
which exact cancellations (insead of resummations) of 
singular terms are expected to occur in physical response 
functions. It appears thus crucial to exploit the underly- 
ing (local-gauge) symmetry and the related Ward identi- 
fied (WI) which implement these exact cancellations, as 
required on physical grounds. In this paper we combine 
the RG approach with the WI to obtain the solution to 
the problem [14]. 

To be more explicit we use the WI to: (i) Establish 
constraints on the RG procedure; (ii) Relate renormal- 



ization parameters to physical quantities; (Hi) Achieve 
the cancellation of singularities in the response functions 
[15]. In this way, the number of marginal and relevant 
running couplings (which are a priori necessary to study 
the IR behavior of interacting bosons for d < 3) is re- 
duced to only one, e.g., the longitudinal two-point ver- 
tex function Tu. In addition, we are able to close the 
equation for Tu, thus providing the exact IR behavior 
for the zero-temperature interacting Bose system. The 
resulting solution is quite different from the Bogoliubov 
one in spite of the coincidence of the linear spectrum. In 
particular, we free the Gavoret and Nozieres results from 
the provisions posed by the occurrence of IR divergences 
[3], and recover the result of Ref. [5] for the anomalous 
self-energy. 

We consider the following action for the Bose system: 
S = J^dr Jdx { V>* (x) \-d r + lA i/>(x) - |(V - iA) ip{x)\ 2 

~\^(x)\*+^(x)X*(x)+r(x)X(x)} , (1) 

where ip{ x ) (with x ~ (r, x)) is a bosonic field obeying 
periodic boundary conditions in the imaginary time t and 
f3 is the inverse temperature (we set h — 1 and m = 1/2). 
In (1) X(x) and {(J,(x), A(x))= A v (x) (y = 0, . . . , d) are 
external "sources" introduced to obtain the correlation 
functions. At the end of the calculation fi(x) recovers the 
constant value fj, of the chemical potential while A and 
A are let to vanish. The interaction potential v is taken 
to be constant in momentum space, as its momentum 
dependence is found to be irrelevant for the IR behavior. 

The action (1) allows for spontaneous broken gauge 
symmetry; in that case, it is convenient to distinguish be- 
tween longitudinal (ipi) and transverse (?/>t) components 
to the broken-symmetry direction, by setting ip(x) = 
ijji(x) + iip t (x) and y^>*(x) = ipi(x) — iipt(x) (with a real 
order parameter). 

By differentiating the free energy F[A v ,Xi] = 
In fT>TjjiT>ipt exp {S}, with i = I, t and A = A; + 
iXt, we obtain as usual the connected correlation func- 
tions, like the "condensate wave function" ipi — 
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{ipi{x)) — SF/SXi(x) and the one-particle Green's func- 
tion Qij — S 2 F/S\iSXj. It is further convenient to 
introduce the Legendre transform of F with respect 
to A, T[A„,i/j io ] = j dxX l (x)ip lo (x) - F[A v ,Xi], whose 
derivatives are the vertex functions ^i 1 ...i n ;v 1 ...v m = 
£("+ m )r / 5tp il0 . . . b^ in0 bA Vl . . . 5A Um associated with the 
one-particle irreducible diagrams of PT. 

In the broken-symmetry phase we keep the value of 
the condensate (ipi(x)) x=0 = tpio fixed. Accordingly, we 
introduce the fields ipi with vanishing averages for van- 
ishing external sources, such that ipi(x) = ipi + il)i(x) 
and 4>t(x) = ipt{x). The mean- field propagators Gij 
are then obtained from the quadratic part of the action 
S Q = -(l/2)Efc^i(-fe)Q«W^(*). wh ere 



H - 



(2) 



with k = (iw„,k) (u> n being a Matsubara frequency). 
In the following we consider the zero-temperature limit 
where to n becomes a continuous variable uj. Enforcing the 
mean- field Bogoliubov condition ipf = /j,/v yields the IR 
behavior Gu ~ (uj 2 + Cgk 2 ) -1 , Git ~ uj(oj 2 + c^k 2 ) -1 , and 
Gu ~ k 2 (w 2 + cjjk 2 ) -1 , where Co = is the mean- 

field value of the sound velocity. This singular IR be- 
havior is associated with the presence of the Goldstonc 
mode. Recall that in the standard ^-representation the 
Bogoliubov propagators (Gn(k) = Q 22 {—k) = Gtt(k) + 
Gll(k)-2iQ lt (k) and Q 12 (k) = Q 21 (k) = -ftt(fc) +Gu{k)) 
share a common (uj 2 + c 2 l k 2 )™ 1 IR behavior, and the nor- 
mal (En (A;) = 2/Lt) and the anomalous (£12 (fc) = ff) 
self-energies satisfy the Hugenholtz-Pines (HP) identity 
£n(0) — £12(0) = 11 [16]. In the ^/^-representation, on 
the other hand, the various propagators have different IR 
behavior since the Goldstone-mode singularity is kept in 
the transverse direction. This choice is crucial to select 
the interaction terms according to their relevance [12]. 

To allow for the RG treatment, it is convenient to 
rewrite the matrix Q(k) in the more general form 



Q(k) 



vu + zuk 2 + uuLU 2 



v it - witw vu + z tt k 2 + u tt uo 



Vl t + WltLO ■ , . 

2 1 „, , ,2 {3) 



where additional terms (running couplings) have been 
introduced with respect to (2). We also introduce run- 
ning couplings for cubic (v ttt , v m , ■ ■ ■) and quartic (v ttt t, 
Vim, Vtttt 1 ■ ■ ■) interaction terms, where the cubic terms 
originate from the presence of the condensate. A pertur- 
bative expansion is then set up, as usual, by regarding 
the quadratic action associated with (3) as the free part 
and the remaining terms as perturbations. In the ab- 
sence of external sources v u and vu vanish by symmetry, 
as shown below. The resulting PT, being massless, is 
plagued by IR divergences already at the one-loop level 
in spatial dimension d < 3 [2,3]. 

Renormalization of the IR divergencies requires a pre- 
liminary power counting for the running couplings. This 



is conveniently done by keeping dimensionless the min- 
imal set of couplings (vu, wu, Ztt) that yields the lin- 
ear part of the Bogoliubov spectrum. We thus rescale 
them and the fields by appropriate powers of cq. In this 
way, [Q tt \ = -2, [Gu] = -1, and [Gu] = 0, where [A] 
stands for the engineering dimensions of A. For sim- 
plicity, from now on we shall omit indicating cq when- 
ever not strictly necessary. We thus have [x] = — 1, 
r] = -1, [ipi(x)] = (d+l)/2, mx)] = (d- l)/2, and 
T t l .u n .!t( k ^ •••)]= ~ni(d+ l)/2 - n t (d- l)/2 + d+ 1. 
The upper critical dimension is d c = 3 [17]. For d < 3 
the running couplings controlling the IR behavior have 
dimensions: [v u ] = [w u ] = [u tt ] = [z tt ] = 0, [v m ] = e/2, 
and [vtttt] — e, with e = 3 — d. Although vu, Vu, and 
Vttt would be strongly relevant, they vanish identically 
for vanishing external sources. 

One could proceed at this point and derive the RG 
equations for the running couplings. As mentioned 
above, however, contrary to critical phenomena in the 
present case of a stable phase, a singular PT has to re- 
sult into finite response functions. It is clear that can- 
cellations of IR divergences in physical quantities signal 
definite connections among the running couplings. In 
the present context these connections stem from the lo- 
cal gauge symmetry and are obtained by examining the 
associated Ward identities [15]. 

In our formalism the WI result from the local gauge 
invariance of the functional T, namely, 



T[A V + d u a(x), R i3 (a(x))ip jo ] = T[A u ,ipi, 



(4) 



Rij(a) being the two-dimensional rotation matrix with 
angle a in the space of the fields ipi and ip t - This equa- 
tion follows from the invariance of the action (1) under 
the gauge transformation tj}{x) — > e la ^ip(x), X(x) — > 
e la(x ^X(x), and A v (x) —> A v (x) + d v a(x) with a(x) real 
function. Taking successive functional derivatives of (4) 
with respect to a, ip i, and A v yields an infinite set of 
WI. For our purposes only the following five WI are rel- 
evant. The first two, which encompass the HP identity, 
are: 

r t i(k)iPio + r t (o) - ik v v hv {-k) = , (5) 

T u {k)^io - r,(0) - ik v T t ; V {-k) = . (6) 

In the limit k v — > they relate the two-point vertices 
to the external sources Aj = T^O), and state that vu = 
rj t (0) and Vu = r«(0) vanish when A, = 0. No gap thus 
appears in the one-particle spectrum. The second couple 
of WI 

r tt /(fci, k 2 )^io + T tt (-k 2 ) - r H (fci + fc 2 ) 

-i(ki) 1/ T t i; V (k2,-k 1 -k 2 ) = 0, (7) 
r t tt(ki,k 2 )4>io - r Jt (-fc 2 ) - T H (ki + k 2 ) 

-i(ki) u TttAk 2 ,-k 1 -k 2 )=0, (8) 
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are the standard WI associated with the continuity equa- 
tion, modified now by the presence of the three-point 
vertices. In the limit k\ = and k 2 — > they yield 
vuttpio = vu and v t tt — 0. The last of our WI is 

^tttt{ki,k 2 ,k 3 )ipi 

-Titt{-k 2 - k 3 ,k 2 ) -Ti tt (ki + k 3 , k 2 ) -T ltt {ki + k 2 , k 3 ) 
-i {k\) v T ttt;v (k 2 , k 3 , -h - k 2 - k 3 ) = , (9) 

from which we obtain Vuttipio/3 = vut for vanishing fc's. 

From the above WI we also relate the running cou- 
plings to the (composite) current vertices and response 
functions. Specifically, in the limit fc^Owe obtain: 



wit = 



1 d 2 T 
4>io dipiodfi ' 



Ztt 



2n s 
(10) 



where n s is the superfluid density (Josephson identity). 
We also have vu — (d 2 T / dipf ) M by its very definition 
[18]. We are left eventually with four running couplings, 
namely, vu, wu, u u , and z tt , whose IR behavior can be 
obtained exactly. 

As a guide to the procedure for obtaining this behavior, 
the RG equations will be evaluated at the one-loop level. 
Exploiting the e-expansion and the minimal subtraction 
technique we get: 



s— 

ds 

dutt 
ds 



cvfi 


dw H 


c Vu Wit 


2^4' 


ds 


' 2^4 


CW H 


s dZtt 


= o, 


2^,4' 


ds 



(11) 



where ,s = k! / k is the scaling factor (k being the normal- 
ization point), c(s) 2 = vii(s)z t t(s) / (vu(s)utt(s) + wit(s) 2 ) 
is the square velocity entering the one-particle propaga- 
tor according to (3), andV>; (s) = ipi K, e / 2 s e / 2 [19]. When 
d = 3 Eqs. (11) are equivalent to those of Ref. [12], pro- 
vided the coupling of Ref. [12] analogous to vutt is iden- 
tified with 3vu/ipf , consistently with the above results. 

Quite generally, the solutions of the coupled equations 
(11) take the form: 

ztt(s) = Ml), wit(s) = — — vii(s), 

Vll(l) 
'D 2 



/ \ n x , Wl t (l) 2 

vu{l) 



W; t (l) 
vu(iy 



This implies that it is sufficient to determine vu(s). 

Although we have derived (12) at the one-loop order, 
we expect them to hold exactly on physical ground owing 
to the identification (10) of the renormalization parame- 
ters with physical quantities. To begin with, the s — > 
value of Ztt is the ratio n s /riQ of finite physical quantities 
(n = ipf being the condensate density) so that diver- 
gences compensate each other in its expression, leading 



to the first of (12). Regarding the second of (12), we 
obtain from (10) for wu and the definition of vu that 
the ratio vu/wu reduces to —2n /(dn /dfi)\ in the limit 
k — > 0. Here (dno/d/j,)\ is the "condensate compress- 
ibility" which has to be finite for thermodynamic sta- 
bility. Finally, from the definition of c(s) 2 and (10) we 
obtain that c(s) 2 reduces to c 2 in the limit s — > 0, where 
now c 2 = 2n s / (dn/d[i)\ is the square of the macroscopic 
sound velocity (n being the density and n = n s at zero 
temperature). By the very stability of the bosonic sys- 
tem, c 2 is free from IR divergences, thus suggesting that 
c(s) is finite and does not scale with s, i.e., c(s) = c = c , 
apart from finite corrections originating from nonsingular 
terms that do not enter the RG flow. Exploiting the first 
and second of (12), we verify that the condition c(s) = 
constant reduces to the third of (12). 

The proof of the above statements is as follows. z u 
has to remain constant by inspection of the WI (6) , which 
shows that the divergence of Zu expected by power count- 
ing is actually not present, since it is related to non di- 
verging quantities. w/ t can instead be identified with vu 
via the WI (7) and (8), which relate vu to Tut and wit 
to r tt; o, respectively, the latter identification being ob- 
tained from the w-derivative of (8). By inspection of the 
leading singular terms to all orders in PT, r tt; o and Tut 
are then found to be proportional to each other [20] ; by 
the same procedure, the invariance of c(s) implied by the 
last of (12) follows from the exact connection between 
the singular parts of r ; oo and Tu , associated respectively 
with u u and vu. 

Determining the IR behavior is thus exactly reduced 
to solving for a single running coupling, for instance vu. 
In particular, at the one-loop order we obtain: 



viljl) 
vu(s) 



1 



vu(l) 



1 + 



2iPi (l) 2 z tt (l) 2 
vu(l)(s-* - 1) 



Ins 



(6 = 0) 

(0 < e < 2) . 



(13) 



2^ (l) 2 ^ t (l) 2 e 

In both cases vu — > as s — > 0, while for e < vu remains 
finite. We show below that the asymptotic behavior (13) 
of vu is actually exact. 

The one-particle Green's function resulting from (12) 
and (13) have the form [19]: 

c n In k 



( 12 ) Gu(k) = 



Gu(k) = 
Gtt(k) = 



8en 2 
dno c 2 



(e = 0) 
(0 < e < 2) 



U) 



1 



dji 4n s uu 2 + c 2 k 2 
1 



c 2 n 



wit 1 _ 

vu z u k 2 + lo 2 /c 2 ' 
1 



2n s 



c 2 k 2 



z tt k 2 +uj 2 /c 2 



where the asymptotic (k — > 0) values of the running cou- 
plings have been identified via (10). Note that the IR 
behavior of Git and Gtt is completely and exactly deter- 
mined by the conditions (12), and is independent from d 
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and the behavior of vu- Instead Gu diverges logarithmi- 
cally as k — > when d = 3 and like k~ e for 1 < d < 3 
[5,6,8]. Accordingly, we find that the anomalous self- 
energy Si2(fc) vanishes like 1/lnfc or like k e , with the HP 
identity now reading £n(0) = /x [5]. In addition to recov- 
ering the result by Gavoret and Nozieres [3] for the prop- 
agators in the ■(/'-representation, W e have also obtained 
the divergent subleading In terms [21]. Such In terms, on 
the other hand, can be shown to disappear altogether in 
the expressions for the density-density, density-current, 
and current-current response functions [5,20]. 

Although vu (s) has been explicitly obtained at the one- 
loop order, we show now that its asymptotic behavior is 
exact. In fact, for d > 3 all interactions are irrelevant, the 
Bogoliubov result is correct, and no longitudinal diver- 
gence appears. For d = 3, all perturbation couplings to 
the Bogoliubov action are marginally irrelevant [12], and 
the result obtained at the one-loop order is stable. For 
1 < d < 3, the one-loop calculation presented above leads 
to a nontrivial fixed point (i.e., Vtttt(s)/n e — > v* ^ 0) 
with vu ~ s e . The WI result v u = v tttt i)i (T) 2 n~ e s e /3 
guaranties that the exponent of the asymptotic behavior 
of vu remains e, irrespective of the actual value v* of the 
nontrivial fixed point, provided it exists. That k e is ac- 
tually the true asymptotic behavior of Tu(k) follows by 
expressing the singular part of Tu in terms of the exact 
T tt i and Q tt : 

T n (k) = B - f - E Stt(q)Gtt(q + k)T ttl (q, -q - k) (14) 
q 

where B coincides with the bare coupling v° t = vu(s = 1) 
apart from finite contributions. The WI (7) allows even- 
tually to express Ttti in terms of Tu itself, thus closing 
Eq. (14) and providing the exact asymptotic behavior 
Tu(k) ~ k e . An analogous procedure was implemented 
in Ref. [5]. 

Some final comments are in order. Contrary to the the- 
ory of critical phenomena, where the WI do not provide 
stringent conditions owing to the genuine divergences of 
physical quantities (such as the compressibility in the 
gas- liquid transition), in our case the relevant physical 
quantities arc finite and the WI imply cancellations in the 
corresponding couplings. In this way, lines of fixed points 
have to appear corresponding to different finite values of 
the thermodynamic derivatives. It is also worth noticing 
that the Bogoliubov mean-field solution for 1 < d < 3 is 
infinitely distant (in the RG sense) from our nontrivial 
fixed point whenever the interparticle interaction is non- 
vanishing. This is true despite the fact that the resulting 
physical picture contains the same low-lying elementary 
excitations. 
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